Covariant gauge on the lattice: a new implementation 
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We derive a new implementation of linear covariant gauges on the lattice, based on a minimizing 
functional that can be interpreted as the Hamiltonian of a spin-glass model in a random external 
magnetic field. We show that our method solves most problems encountered in earlier implementa- 
tions, mostly related to the no-go condition formulated by L. Giusti, Nucl. Phys. B 498, 331 (1997). 
We carry out tests in the SU(2) case in four space-time dimensions. We also present preliminary 
results for the transverse gluon propagator at different values of the gauge parameter £. 
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I. INTRODUCTION 

The behavior of Green's functions in the infrared (IR) 
limit of Yang-Mills theories is of fundamental importance 
for the understanding of the low-energy properties of 
quantum chromodynamics (QCD), in particular for the 
problem of quark and gluon confinement Since the 
evaluation of these functions depends on the gauge condi- 
tion, it is important to consider different gauges in order 
to obtain a clear (possibly gauge-independent) picture 
of color confinement. Needless to say, this investigation 
should be at the nonperturbative level. 

A nonperturbative study of the QCD propagators and 
vertices from first principles is possible using lattice sim- 
ulations. Of course, on the lattice, the finite size of the 
system corresponds to an IR cutoff ~ 27r/L, where L 
is the lattice size. Thus, a numerical study of Green's 
functions in the IR limit usually requires a careful ex- 
trapolation of the data to the infinite- volume limit Q. 
Another possible limitation for the simulations is the dif- 
ficulty in finding an efficient numerical implementation 
of a given gauge condition. For this reason, most numer- 
ical studies of Green's functions have been restricted to: 
Landau gauge [|[, Coulomb gauge [U, A-gauge (a gauge 
that interpolates between Landau and Coulomb) [5( and 
maximally Abelian gauge 6]. On the other hand, among 
the various gauge conditions that are very popular in 
continuum studies, the so-called linear covariant gauge 
— which is a generalization of Landau gau ge — pr oved 
quite hostile to the lattice approach 0, d, HHH EH EH • 

Let us recall that, in the continuum, the linear covari- 
ant gauge is defined by 



d M A b Jx) = A»(x) , 



(1) 



where A^(x) is the gluon field and the real- valued func- 
tions A b (x) are generated using a Gaussian distribution 

P[A b (x)] ~ exp{- £ [A b (z)] 2 /(2£)} ■ (2) 



The Feynman gauge corresponds to the value £ = 1, while 
the Landau gauge is obtained in the limit £ — > 0. 

In the Landau case, the gauge condition d^A^x) = 
can be obtained by minimizing the functional 



£ LG {A*} ex J d^x ]T [(A%(x)} 



(3) 



with respect to the gauge transformations {g(x)}. In 
Ref. 0] it was shown that a similar minimizing functional 
£lcg{A 9 } for the linear covariant gauge does not exist. 

The no-go theorem proven in [7[ can of course be 
avoided by relaxing its hypotheses. The first possibil- 
ity, explored in that reference, is to consider a different 
gauge, i.e. Fid^A^x) - A b {x)} = with F[0] = 0, for 
which a minimizing functional exists. In particular it was 
shown that minimizing the functional 



d 4 s 



(j, . b 



{[d ll A b fl (x)-A b {x)]*} 



implies the stationarity condition 

Dl b d v [d^{x)-A\x)} 



0. 



(4) 



(5) 



where Z?° fc is the covariant derivative. 

As noted previously 0], this method presents several 
problems. First of all, one can introduce spurious solu- 
tions, corresponding to F[s] = for s ^ 0. In the above 
case, these solutions are the zeros of the operator D^d^ . 
Also, the second derivative of this functional does not cor- 
respond to the Faddeev-Popov operator M. = —d tl D^ b 
of the usual linear covariant gauge. Finally, the lattice 
discretization of the above functional is not linear in the 
gauge transformation {g(x)}. This makes the numerical 
minimization difficult and one has to rely on a specific 
discretization of the minimizing functional [|| in order to 
make the lattice approach feasible. 

A second possibility, recently presented in [ll[, is 
based on avoiding the use of a minimizing functional 



2 



£lcg{A 9 }, i.e. on considering a lattice definition of the 
linear covariant gauge that coincides with the perturba- 
tive definition in the continuum. In this case, one first 
fixes the gluon field to Landau gauge, i.e. the transformed 
gauge field satisfies d^A* Jx) = 0. Then, one solves the 
equation {d^D^cj) ) (x) = A b (x) and uses 4> c ( x ) as a 
generator of a second gauge transformation. For small 
4> c (x), one then has that the gauge-transformed gluon 
field A' Ax) satisfies the condition 



is the minimizing functional for the lattice Landau gauge. 
Here, the link variables U^x) and the site variable g(x) 
are matrices belonging to the SU(iV c ) group (in the fun- 
damental representation). We also indicate with the 
real part of a complex number and with Tr the trace in 
color space. Indeed, the lattice linear covariant gauge 
condition can be obtained by minimizing the functional 



£ 



LCG 



{U 9 ,g} = £ LG {U 9 } + $Tr}2 ig(x)A(x) . (9) 



d,A'l(x) = d, (4 + D b ^ c ) (x) = A b (x) . (6) 

The main problem in this case is that the method is 
correct only for infinitesimal gauge transformations, but 
usually 4> c (x) is not small in a numerical simulation. As 
a consequence, one finds dH that the distribution of 
d^A'^x) does not agree completely with the Gaussian 
distribution of A b (x). Moreover, in linear covariant gauge 
the longitudinal gluon propagator Di(p 2 ) should satisfy 
the relation p 2 Di(p 2 ) = £. With this approach one finds 
that, for small momenta, this is not the case [111 ]. 

Here we present a new implementation of the linear 
covariant gauge on the lattice that solves the problems 
illustrated above, afflicting earlier methods. The paper 
is organized as follows. In the next section, we describe 
our new implementation and we report some tests of the 
algorithm used for the numerical gauge fixing. In par- 
ticular, we show that Di(p 2 ) is well described by £/p 2 , 
where £ is the gauge parameter, and discuss discretiza- 
tion effects. We also present preliminary results for the 
momentum-space transverse gluon propagator D t (p 2 ) for 
different values of £. Finally, we present our conclusions. 



II. A NEW IMPLEMENTATION 

In order to find a new implementation for the linear 
covariant gauge on the lattice we notice that, when min- 
imizing a functional £{A 9 }, the gauge condition is given 
by the first variation of £, i.e. 



ic d£ dA 
S£= dA^ S9 







(7) 



Ref. [?J has proven that there is no functional £ {A 9 } lead- 
ing to Eq. (fTJ). Nevertheless, one can remove an implicit 
hypothesis of the no-go condition, i.e. that the gauge 
transformation {g(x)} appears in the minimizing func- 
tional in the "canonical" way A 9 . Thus, we may look for 
a minimizing functional of the type £lcg{A 9 , g} instead 
of simply £ LC g{A 9 }. 

If one recalls that solving the system of equations 
Bip = £ is equivalent to minimizing the quadratic form 
rpBip/2 — V'Ci then it is clear that we should have 
£ L cG{A 9 ,g} ~ £ LG {A 9 } - gA, where 



£ LG {U g }=-^TrJ2 9{x)U ll {x)g\o 



To prove that this is the right functional we can consider 
a one-parameter subgroup g(x,r) = exp [iTj b (x) \ b ~\ 
of the gauge transformation {g(x)}. Here, A & are the 
traceless Hermitian generators of the Lie algebra of the 
SU(_/V C ) gauge group. They also satisfy the usual normal- 
ization condition Tr (A b A c ) = 2S bc . Then, it is easy to 
check that the stationarity condition d T £LCG{ T = 0) = 
[for all 7 b (a;)] implies the lattice linear covariant gauge 
condition 



V-A b (x) =J2 A U X ) 



K(x- 



- e„) = A b (x) . (10) 



(8) 



Here, A b (x) = Tr A(x) \ b and, similarly, A^x) = 
Tr A^(x) X b with 

A»(x) = (2i)-i [U»(x)-Ut(x)] traceless . (11) 

Note that Eq. dTOj) implies that Y, x A 6 (a;) = 0. 

At the same time, one can verify that the second term 
on the r.h.s. of Eq. © does not contribute to the sec- 
ond variation of the functional £lcg{U 9 , g} with respect 
to the parameter r. Thus, this second variation defines 
a matrix M. that is a discretized version of the usual 
Faddeev-Popov operator —d-D. Let us recall that, due to 
the gauge condition ([T]) above, in linear covariant gauge 
one has in general —d ■ D ^ —D ■ d. 

Finally, let us note that the functional £lcg{U 9 ', g} 
in Eq. ([9]) is indeed linear in the gauge transforma- 
tion {g(x)}. Also, one can interpret the Landau-gauge 
functional £lg{U 9 } [see Eq. ([5]) above] as a spin-glass 
Hamiltonian for the spin variables g{x) with a random 
interaction given by U^(x). Then, our new functional 
£lcg{U 9 , g} corresponds to the same spin-glass Hamil- 
tonian when a random external magnetic field A(x) is 
applied. 

We have performed some numerical tests with the func- 
tional ©, using the so-called stochastic- overrelaxation 
algorithm [3]. To this end we considered the 4d SU(2) 
case at /3 = 4, for V = 8 4 and 16 4 , with £ = 0.001,0.1 
and 0.5. The numerical gauge fixing seems to work very 
well in these cases (see plot in Figure [T]). 

We also checked that the quantity p 2 Di(p 2 ) is constant 
within statistical fluctuations in all cases considered. For 
V = 16 4 and £ = 0.5 a fit of the type a/p b for Di{p 2 ) 
gives a = 0.502(5) and b = 2.01(1) with a x 2 /dof = 1.2. 
Similar fits have been obtained in the other cases. 
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FIG. 1: Convergence of the numerical gauge fixing. Here we 
report the value of A = 6 [V- A b (x) — A h (x)] 2 as a function 
of the number of iterations n for (3 = 4, £ = 0.5 and V — 8 4 . 
Note the logarithmic scale on the y axis. 

A. Discretization effects 

The minimizing functional ([9|) in principle solves the 
problem of fixing the gauge condition (jTUJ) . From the nu- 
merical point of view, however, one has to recall that the 
gluon field is bounded, at least when using the standard 
(compact) discretization. Thus, while our method works 
wery well when the functions A b (x) are generated using 
a bounded distribution, care must be taken in the usual 
implementation of the linear covariant gauge, where the 
functions A b (x) [see Eqs. (JTJ and (jSJ)] satisfy a Gaus- 
sian distribution, i.e. they are unbounded. This can give 
rise to convergence problems [f3 ] when a numerical im- 
plementation is attempted. Of course the problem gets 
more severe when £ is larger. This is a common problem 
of all lattice realizations of the linear covariant gauge. 

Actually, in order to obtain the correct continuum limit 
the functions A b (x) are generated on the lattice, 
in the SXJ(N C ) case, from a Gaussian distribution with 
width y/a = y/2N c £/f3, instead of the width V?- Note 
that, for (3 = 4 and N c = 2, as in the previous section, 
one has a = $. Thus, for a given value of £ one can obtain 
a sufficiently small value for a by considering large values 
of the lattice coupling (3. However, if j3 is too large the 
physical volume will be too small (for a given lattice size) 
and one cannot really probe the IR limit of the theory. 
On the other hand, for (3 < 2N C the lattice width y/a is 
even larger than the continuum width y/£. Note that the 
situation is probably better in the ST/ (3) case, since one 
has cr = £ for /3 = 6, which corresponds to a reasonably 
large value of the lattice spacing a. 

Of course, one can try to use different discretizations 
of the gluon fields in order to improve the convergence of 
the minimizing algorithms. Besides the usual discretiza- 
tion (fTTj) . we also did some tests with the "angle" projec- 
tion [1 51 ] and with the recently introduced stereographic 
projection (or modified lattice Landau gauge) [16] , us- 
ing the so-called Cornell method [l3j]. In the last case 
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FIG. 2: Transverse gluon propagator D t (p 2 ) as a function 
of the momentum p (both in physical units) for the lattice 
volume V = 16 4 , f3 = 2.3 and £ = (+), 0.05 (x ), 0.1 (*). 

the gluon field is in principle unbounded even for a finite 
lattice spacing. Indeed, our tests show that with this 
discretization one is usually able to simulate at slightly 
larger values of £, for a given lattice volume V and lattice 
coupling (3, compared to the other two cases. In partic- 
ular, we tested these three discretizations using V = 8 4 , 
£ = 0.01, 0.05, 0.1, 0.5, 1.0 and (3 = 2.2, 2.3, . . . , 2.9, 3.0. 
We found that, while the usual discretization and the 
"angle" projection have problems with £ > 0.5 already 
at (3 — 2.9, the stereographic projection allows one to 
simulate with £ = 1 for (3 > 2.5. Note that £ = 1 and 
j3 = 2.5 corresponds to a = 1.6. 



B. Transverse gluon propagator 

Using the minimizing functional shown above and the 
stereographic projection, we have simulated at (3 = 2.2 
and j3 = 2.3 for the lattice volumes V = 8 4 , 16 4 and 
24 4 , for several values of the gauge parameter £ in the 
SU(2) case. As a test of the gauge-fixing method, we 
have checked that the quantity Di(p 2 )p 2 /a, which should 
be equal to 1, has a value of 0.999(2) when averaged over 
all data D/(p 2 ) produced. Preliminary results for the 
transverse gluon propagator D t (p 2 ) as a function of the 
momentum p are shown in the Figures [2] and [3] There is a 
clear tendency of getting a more suppressed IR propaga- 
tor when the lattice volumes increases, as in the Landau 
case, and also when the value of £ increases. The latter 
result is in agreement with Ref. [1J|- Also, the extrapola- 
tion to infinite volume, for a given (3 and a fixed value of 
£, seems in this case even harder than in Landau gauge. 
Indeed, as V — > oo, the number of sites characterized by 
a large value for the function A b (x) increases, making the 
convergence of the gauge-fixing method more difficult. 

One should recall here that, at the perturbative level, 
the gluon field A and the gauge parameter £ are (mul- 
tiplicatively) renormalized by the same factor Z3, i.e. 

1 /2 

Ab = Z 3 Ar and £b = Z^^r, where as usual B 
and R indicate bare and renormalized quantities respec- 
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FIG. 3: Transverse gluon propagator D t (p 2 ) as a function 
of the momentum p (both in physical units) for the gauge 
coupling £ = 0.05, j3 = 2.3, with the lattice volumes V = 
8 4 (+),16 4 (x) and 24 4 (*). 

tively. On the lattice this implies that, in the scaling 
region, data obtained for two different values of /3, e.g. 
(3\ and /?2, can be compared and should give the same 
(renormalized) propagator only if the multiplicative fac- 
tor Rz = Zs(Pi) I Z 3 (f3 2 ) relating the propagators (see 
pjj for the case of Landau gauge) also relates the gauge 
parameters £i and £2- Since the value of Rz is not known 
a priori, one has to use some type of "matching proce- 
dure" in order to find pairs of parameters (/?, £) yielding 
the same continuum renormalized propagators. This is 
left to a future work [l8| . 

III. CONCLUSIONS 

We presented a new lattice implementation of the lin- 
ear covariant gauge by using a minimizing functional 



£lcg{U 9 ', g}. Tests have been done for the SU(2) case in 
four space-time dimensions. This approach solves most 
problems encountered in earlier implementations and en- 
sures a good quality for the gauge fixing with a ratio 
Di(p 2 )p 2 1 X ~ 1 for all cases considered. We have also re- 
ported preliminary results for the transverse gluon prop- 
agator D t (p 2 ). The only open problem is how to extend 
these simulations to large lattice volumes, in order to 
probe the IR limit of the theory when the gauge param- 
eter £ is also large. We stress that infinite-volume results 
for D t {p 2 ) would of course be v ery important for compar- 
ison with analytic studies [pj [2fj | in the continuum. In 
any case, as mentioned above, the discretization effects 
are probably less severe for the SU(3) group compared 
to the SU(2) case. We are currently simulating other 
values of f3 and £ in the 4d SU(2) case and considering 
simulations also of the SU(3) group and of the 3d case 

Finally, having a minimizing functional for the linear 
covariant gauge, which extends in a natural way the Lan- 
dau case while preserving all the properties of the con- 
tinuum formulation, allows a numerical investigation of 
the effects of Gribov copies on Green's functions for the 
case £ 7^ 0. These studies will make possible a compari- 
son to recent analytic results [2(| obtained in the limit of 
small gauge parameter £. They could also be important 
for understanding how the so-called (Landau) Gribov- 
Zwanwiger confinement scenario Q should be modified 
in the general covariant gauge. To this end, results ob- 
tained for small values of £, which can be easily obtained 
using our new implementation, could already be relevant. 
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